Kondo resonant spectra in coupled quantum dots 
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The Kondo effect in coupled quantum dots is investigated from the viewpoint of transmission 
spectroscopy using the slave-boson formalism of the Anderson model. The antiferromagnetic spin- 
spin coupling J between the dots is taken into account. Conductance G through the dots connected 
in a series is characterized by the competition between the dot-dot tunneling coupling Vc and the 
level broadening A in the dots (dot-lead coupling). When Vb/A < 1, the Kondo resonance is 
formed between each dot and lead, which is replaced by a spin-singlet state in the dots at low gate 
voltages. The gate voltage dependence of G has a sharp peak of 2e 2 /h in height in the crossover 
region between the Kondo and spin-singlet states. The sharp peak of G survives when the energy 
levels are different between the dots. When Vc/A > 1, the "molecular levels" between the Kondo 
D [ resonant states appear; the Kondo resonant peaks are located below and above the Fermi level in 

the leads at low gate voltages. The gate voltage dependence of G has a broad peak, which is robust 
against J. The broad peak splits into two peaks when the energy levels are different, reflecting the 
formation of the asymmetric molecular levels between the Kondo resonant states. 
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O ; I. INTRODUCTION 

Recently the. Kondo effect has been observed in semiconductor quantum dots connected to external leads by tunnel- 
ing barriers.Ercl T-he-Kondo effect makes a resonant state at the Fermi level in the leads when the number of electrons 
in the dot is oddJ3il3 The resonant width is given by the Kondo temperature Tk . This results in (i) the unitary limit 
of the linear conductance through the dot G = 2e 2 /h at low temperatures (T <C Tk) (Ref. ||) and (ii) the zero-bias 
peak of the differential conductance dI/dV s d with the width of Tk under finite source-drain voltagesfl 

In this paper, we theoretically examine the Kondo resonant state in double quantum dots connected in a spjj 6 ^ as 
' will be shown in Fig. 1. We consider one level in each quantum dot, often referred to as an "artificial -atom" £311-3 At 
T 3> Tk, the transport properties of such systems have been investigated by several experiments,E3~E3 from a pfiWi 
point of the formation of an "artificial molecule." The molecular orbitals between the two dots have been observed jUrEij 



which reflect the coherent coupling between the dots (Vc in Fig. 1). When the doub 



it system accommodates two 



electrons, the coherent coupling and intra-dot interaction U make many-body statesEj'El The ground state is a spin 
singlet. The excitation energy to a spin-triplet state is given by 



V(U/V C ) 2 + 16 - (U/Vc) (1) 



when U ^ Vc- The ground state is an "entangled" state between two localized spins, which has attracted new interest 
for the application to quantum computing .EU At T <C Tk, the two elements Vc and J would lead to rich structures 
of the Kondo resonant state in this system. 



A-,useful tool to elucidate the Kondo resonance is the slave-boson formalism of the Anderson model with U 



2 J~E£l In the case of a single quantum dot, the dot state is spin- up /||0), spin-down /||0), or empty b^\0) (/|/-| + 
fffl + = !)• I n the mean-field approximation, b and are replaced by a real c-number, b^ = b = &o- The density 
of states for pseudo fermions per spin is0 

. . Im 1 1 A 

P f{uj)= » =r = ^ — . 2 

TT LU - E + iA 7T (w - E) 2 + A 2 



1 



In the Kondo regime where one electron exists in the dot (6q <C 1, /j/j + flfi — 1), the density of states /0/(w) is 

half filled and, hence, the center of the resonant state E is matched with the Fermi level. The width A is the Kondo 
temperature 

T Kil = - exp(7r£7A), (3) 

where A is the level broadening, 7rp|V| 2 , with the density of states p in the external leads and the dot lead tunneling 
coupling V. E*jjs the renormalized dot-level, E* = V g + A/ir In (Dtv/A) with the gate voltage V g and the bandwidth 
D in the leads.c3 Therefore this method is suitable to examine how the Kondo resonant state is formed. When the 
number of electrons in the dot is less than unity (valence fluctuating regime), the center of the resonance moves 
above the Fermi level. The resonant width increases from Tj( to A, as the number of electrons in the dot decreases. 
However, the mean-field approximation is not quantitatively accurate in the valence fluctuating regime. 

The KcmdQ effect on the transport, properties in coupled dot systems have been studied by several theoretical 
methods. E3~ a In our previous paper ,E3 we have applied the slave boson mean-field theory to the double dot system in 
the absence of the antiferromagnetic spin-spin coupling J. We have shown that the Kondo resonance is determined 
by the competition between the dot-dot tunneling coupling Vc and dot-lead coupling A. The resonant state has two 
peaks below and above the Fermi level for Vc/A > 1, whereas it has a single peak at the Fermi level for Vc/A < 1. 
In consequence, the linear-conductance G through the double dots as a function of the gate voltage is qualitatively 
different between the two cases. „ 

Extending our calculations, Georges and Meirtj have pointed out the importance of the spin-spin coupling J. At 
low-gate voltages, the coupling becomes relevant and makes a spin-singlet state in the two dots. Then the Kondo 
resonance between dots and leads is destroyed. Izumida and SakaiEj have examined the same problem using the 
numerical renormalization group (NRG) method. A-similar situation was investigated for dilute magnetic systems, 
which is known as the two impurity Kondo problem. E3e3 When the antiferromagnetic coupling J, the origin of which 
is the Ruderman-Kittel-Kasuya-Yosida interaction, increases from zero, electronic states of two magnetic impurities 
exhibit the first-order phase transition from the Kondo state to the spin-singlet state at J ~ Tk,\- The m agn etic 
susceptibility diverges at the transition. The phase transition changes into a crossover in the case of Vc ^ QJ^H In 
coupled quantum dots, J is given by-En. (|l]) whereas Tk,i increases as the gate voltage increases as in Eq. (||). 

As indicated by the previous workL3LJ, the interplay of A, Vc, and J results in various transport properties. They 
are observable by experiments since these parameters can be controlled in quantum dot systems. The physical origins 
of the transport properties, however, are difficult to understand by the NRG methodEJ and by the conventional 
argument using the phase shift.ESI In this paper, we look directly at the Kondo resonant state in the presence of 
spin-spin coupling J. The effect of J can be understood in terms of the transmission spectrum, just as in our previous 
work in the absence of J. We will show that even for J ^ 0, electron transport is characterized by the competition 
between the dot-dot tunneling coupling Vc and the dot-lead coupling A. When Vc/A < 1, the coupling J is essentially 
important. G has a sharp peak of 2e 2 /h in height, as a function of gate voltage, which stems from the coexistence 
of the Kondo coupling and spin-singlet coupling. When Vc/A > I, the transport properties are robust against J. 
These results are explained by only one criterion, whether the transmission spectrum has double peaks. Hence our 
argument is much simpler than the one based on the phase shiftO In addition, the transmission spectrum by itself 
is directly_observable in the differential conductance under finite source-drain voltages, as suggested by Aguado and 
Langretho Further, we generalize our calculations to asymmetric double dot systems. Some preliminary results have 
been reported in Ref. ^4|. 

The organization of this paper is as follows. In Sec. [n], t ne m odel Hamiltonian for a coupled quantum dot system 
is introduced. We explain the slave boson mean-field theory and derive the expression of the conductance through 



the dots. In Sec. Ill, we briefly review electron transport through symmetric double dots in the absence of J. The 



transmission spectra are discussed in detail. In Sec. [TV| , we investigate electron transport through coupled dots in 
the presence of J. In Sec. |V|, the electron transport through asymmetric double dots is examined. Conclusions and 
discussion follow in Sec. 



II. MODEL AND METHOD 



A. Model 

Let us consider coupled quantum dots as shown in Fig. |l|. Two dots couple to each other with Vc, and to external 
leads with V. Each dot has a single energy level E a (a — L,R). We denote the energy difference, El — Er, by AE. 
A common gate V g is attached to the dots to control E a ; we define V g = (El + Er)/2. We assume that the intradot 
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Coulomb interaction U is sufficiently large so that (i) the double occupancy of electrons in each dot is forbidden, but 
(ii) the antiferromagnetic spin coupling exists between the quantum dots due to the virtual double occupancy in a 
dot, JS L ■ S R with J (= AV§/U) [see Eq. (g)] and the spin operator S a (a — L, R). The interdot Coulomb interaction 
is neglected. The temperature is set to be zero. 

We adopt the (N = 2)-fold degenerata-Anderson model with the antiferromagnetic SU(JV) spin interactions, rep- 
resented by the slave-boson formalism Bj^LJ We introduce the slave boson operator b^ a , which creates an empty state 
and fermion operator /^ m , which creates a singly occupied state with spin m = —j, —(J — 1), ...,j — l,j(2j + 1 = AT) 
in dot a — L,R, under the constraint of Q a = n a ^k bLb a = 1 with n a = ^2 m famfam- The annihilation operator of 
an electron in dot a is rewritten as C am = &^/ am £3~E3 

The Hamiltonian is 

H = Ho + Hdot — lead + ^dot-dot + Uj + E A «(<3«-i)- (4) 

a=L,R 

The first term on the right-hand side of Eq. (Q) represents electrons in the leads and dots; 

Ho = E E ( k )^km^k m + E E a p m J am . (5) 

a = L,R a = L ; R 

k, m 7n 

The operator c^ fem creates an electron in lead a with energy E(k) and spin m. The second and third terms on the 
right hand side of Eq. (|J) describe the dot-lead tunneling coupling and the dot-dot tunneling coupling, respectively 

Hdot-lcad = -7= E C lkmf^ h i + H.C., (6) 
a,k,m 

^dot-dot = ^ E b ifLmfL n b R + H.C. (7) 

m 

The term Hj gives the antiferromagnetic spin-spin interaction 

Hj = Jf E fRmhrnfLfRn + H.C. (8) 
m.n 

The constraint of Q a — 1 is taken into account by the last term on the right-hand side of Eq. (||) with the Lagrange 
multiplier A Q . 



B. Mean-Field Approximation 

Let us now make the following two assumption to treat the Hamiltonian (|f|) within the mean- field theoryiaB First, 
the slave-boson operators b a and b' a are replaced by a constant c- number b a . This approximation is psjaet in the limit 
of N — ► oo when J = 0; it corresponds to the calculation to the lowest order in the 1/N expansionEJeJ Second, we 
also treat the spin-spin coupling (|8j) by the mean-field theory. It is decoupled by a parameter k, 

m 

which characterizes the antiferromagnetic order between the spins in the dots. Then we obtain the following mean-field 
Hamiltonian; 

' E L V c +n\( f Lri 
^ W^afcrn^fcrn -r / ,\j Lm j Rm ) y ~ 

a = L,R 

k. m 



n =j: E { k)ci km c akm + E(/L/D ~*l Y ){h 

=l,r m \ C» R / \ 

, rn 

■-7ft E 6 «(4m/«m + /L^m)+ E A « ^ ~ ^ + J ^ > ^ 

™ ^ ™ a— L,R 
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where E a = E a + X a and Vc = b~i}>RVc with b a = b a /^/~N. 

From this Hamiltonian, we can calculate the ground state energy. Minimizing the ground state energy, we determine 
five parameters b a , X a (a — L, R) and n self consistently. (See Appendix A.) In general, the self-consistent equations 
have two solutions; one is A Q = (b a — 0) and k = J/2, and the other is A a ^ 0. The former means that electrons in 
the dots are isolated from the leads and make a spin-singlet state. The latter corresponds to the coherent couplings 
between dots and leads. We determine the ground state by comparing the energies of these solutions. 



C. Conductance 



Next we derive the expression of the conductance G. Current / through the coupled dots under a source-drain 
voltage V s d is given by a charge transfer through the central tunneling barrier 



-ieVc 



E(<4 



L m CR m ) 



(4 



Rm C Lr 



(11) 



where (...) represents the Keldysh-Green function. By the mean-field approximations, (c Lm CR m ) = 6L^(/^ m /ij rra ), 
where (f Lm fR m ) can be calculated using the Hamiltonian (10). (See Appendix B.) 

Evaluating the Keldysh-Green functions, we obtain the expression of the conductance G = limy ad ^o 1 7 Vsd> 



2e 2 

G = — 7> = 0), 



(12) 



formally in the same form as the Landauer formulae£3 The transmission probability T{uj) through the two dots is 

4A L A^ 2 



2» 



(u-E L - tA L )(uj - E R - iA R ) -t 2 (uj-E l + tA L )(uj - E R + iA R ) - t 2 



(13) 



with A Q = 6 2 A (a = L, R) and t — Vc + n. We choose u> = at the Fermi level in the leads. 



III. SYMMETRIC DOUBLE DOTS WITHOUT J 



We begin with a symmetric double dot system El = Er in the absence of the spin-spin coupling J — 0. Because 
of the symmetry between the two dots Bl = t>R = b,riL = tlr = n/2, Xl = Xr, and El — Er = E. k — due to the 
absence of J. Then the expression of the transmission probability T(ui) in Eq. (fL3ft is simplified to 



T(w) 



4A 2 V 2 



(oj-iE + Vc)) +A 2 



u-(E-Vc)) +A 2 



(14) 



with A = b 2 A and Vc = b 2 Vc- The parameter Vc characterizes the effective dot-dot coupling. Note that Vc/A = 
Vc/A. ' 

The line shape of T{uj) is determined by the parameter of Vc/A£3 In Fig. 0, T(u>) is plotted as a function oiu>/A. 
When Vc/A < 1, T{u>) has a single peak at u> = E [Fig. [|(a)], whereas when Vc/A > 1, T(u>) has two maxima at 

lu = E T \fV~c ~ ^ 2 [ Fi §- i( c )l- When V c / A = 1, T(lu) has a flat-topped single peak as shown in Fig. |(b). 

The dotted line in Fig. ||(a) represents the conductance as a function of the gate voltage when Vc/A < 1. For 
sufficiently low gate voltages, the Fermi level in the leads is located at the peak of T(ui) [the case of V g /A = —2 is 
indicated by the arrow in Fig. ||(a)]. In this case, each dot accommodates one electron and forms the Kondo resonant 
state with conduction electrons in a lead. The electron transport is determined by the hopping between the two Kondo 
resonant states. Thus G is proportional to Vq and independent of the gate voltage. As the gate voltage increases, the 
position of the Fermi level shifts downwards [the cases of V g /A = and 1 in Fig. ||(a)]. The number of electrons in 
the dots n decreases monotonically from 2 to (valence fluctuating regime). As a result, G decreases with increasing 
the gate voltage. 

The dotted line in Fig. [|(a) represents the conductance as a function of the gate voltage when Vc/A > 1. For low 
gate voltages, the Fermi level in the leads is located at the center of T(u) [the case of V g / A = —1 in Fig. 0(c)], and 



4 



thus the conductance G is considerably suppressed. As the gate voltage increases, the Fermi level shifts downwards 
[the case of V g /A = —4 in Fig. @(c)], which increases G. At a gate voltage of V g « 0, the Fermi level is just on the 
bonding peak and, thus, the conductance G has the maximum value of 2e 2 /h. On increasing the gate voltage further, 
the Fermi level shifts downwards further [the case of V g /A = 5 in Fig. 0(c)], which leads to a decrease in G. 

In summary, when Vc / A < 1, the electron transport is characterized by electronic states in each dot below the 
Kondo temperature whereas, when VpJ A > 1, it is characterized by the formation of molecular levels of the Kondo 
resonant states. In the recent workjlj Aguado and Langreth have argued that these molecular levels are directly 
observable in the differential conductance under finite source-drain voltages. 

We next mention V g dependence of the peak width A of the Kondo resonances. When Vc/A <C 1, 

?* <y, < o) fl51 

A VF = A/2 (| + f arctan( VA)) (0 < V g ) ' [ °> 

where Tk is the Kondo temperature for the coupled dots, 

V c V c 



T K = T K<1 exp arctan j /y/l + (V c /A) 2 . (16) 

Tk,i is given by Eq. (||). Tk increases exponentially with V g or Vc-S When Vc/A^> 1, 

f Tk (V g < -Vc/2) 

A = I 4^(V g + Vc/2) (-Vc/2 <V g < Vc/2) . (17) 

I Avf (Vc/2 < V g ) 

(See Appendix C.) 

IV. SYMMETRIC DOUBLE DOTS WITH J 

Now we discuss a symmetric double dot system, El = Er, in the presence of the spin-spin coupling J. When 
the coupling J is strong enough, the Kondo coupling between a dot and lead is destroyed and a spin-singlet state 
appears between the dots. First we make a rough estimation of J = J e , w here the crossover between the Kondo and 



spin-singlet states takes places. The ground-state energy e gs in Eq. (All) is given by 



2~ k 2 2EA , . 

e gs = ^--A- y -^, (18) 

where Al,j? = A, = E, and til^r = n/2. Here we have used the self-consistent equations in Appendix A. The 

last term is negligible when n ~ 2. We consider two simple situations; (i) Kondo state without spin-singlet ordering 
(k = 0), e gs = V g — -A, and (ii) spin-singlet state without the Kondo effect (A = and n = J/2), e gs — V g — J/4. 
The crossover between (i) and (ii) takes place at J c /A = 8/tt = 2.54. If Vc — > 0, A — > T/^j_and thus J c /Tk,i — 2.54. 
This is consistent with the NRG result of the two impurity Kondo problem, J c /Tka = 2.20 When Vc/ A = 0.11 and 
J c /A = 2.1.x 10~ 3 (£//A = 21.2), J c /A = 2.55 is in good agreement with the recent NRG calculation of the coupled 
dot systemEJ that gives J c /Tk,\ = 2.65. 

In coupled quantum dots, J is given by AVq/U. U/A = 4 x 10 2 is assumed in this paper. We investigate the gate 
voltage dependence of G with fixed J. With decreasing V g , the Kondo effect becomes weaker, as indicated in Eq. (||), 
and thus J is relatively stronger compared with A. At sufficiently low gate voltages, a spin-singlet state appears that 
destroys the Kondo effect. 

A. Transmission probability 



In the present case, the transmission probability T(uj) is given by 

4A 2 * 2 



(co -(£ + <))' + A 2 (co- (£-<)) 2 + A 2 



(19) 



where t = Vc + k, Vc = b 2 Vc, and A = 6 2 A. This result is the same as the one without J [Eq. (jlj)] if the dot-dot 
coupling Vc is replaced by t. Then the line shape of T(lo) qualitatively depends on whether t/A = Vc/A + k/A is 
larger than unity. 
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B. Weak dot-dot coupling case 



First we discuss electron transport when Vc/A < 1. The solid line in Fig. [|(a) represents V g dependence of G, 
while the solid and dotted lines in Fig. |^(b) represent V g dependence of k and A, respectively (Vc/A = 0.3, J/A = 
9.0 x 10~ 4 ). For sufficiently low gate voltages, k — J/2 and A = 0, indicating a complete spin-singlet state, which 
results in G = 0. With increasing V g , the Kondo effect is stronger. As a consequence, A becomes larger, whereas k 
becomes smaller. At high gate voltages, G is identical to the one without J (dotted line). Here the spin-spin coupling 
J is weak enough compared with the Kondo effect. Between these two regions, G has a sharp peak of 2e 2 /h in height. 
At this peak, both of the Kondo effect A and spin-sincdet coupling k coexist [see the inset in Fig. ||](b)]. This peak 
structure of G is consistent with the NRG calculations.^ 

The calculated results of the G-V g curve can be understood by Eq. ([L9|). At high gate voltages, k/A < 1 so that 
t/A = Vc/A + k/A ~ Vc/A; the plot of G is similar to the one of J = 0. With decreasing V g , k/A increases and 
thus t/A becomes larger. On the right hand side of the sharp peak of G, t/A < 1 and hence T(oj) has a single peak. 
At k/A = 1 — Vc/A, t/A = 1 and T(oj = 0) = 1, resulting in G = 2e 2 /h. With decreasing V g further, split peaks 
of T(ui) suppress G. At a certain gate voltage, the Kondo states disappear and a spin-singlet state appears (n = J/2 
and A = 0). 

From the results given above, the origin of the peak of G is explained as follows. The Kondo resonant state is 
formed by the coherent coupling between the dot and the lead while the spin-singlet state is formed by the coherent 
coupling between the two dots. When these two states coexist, the left and right leads are connected via a coherent 
channel. Electrons transfer through this channel and as a result, T(io = 0) = 1. 

We find that the width of the sharp peak is roughly given by Vc- This implies that the crossover region between 
the spin^ainglet state and the Kondo state becomes wider as Vc increases. It is consistent with the results by Sakai 
et a/.oEl for the two-impurity Kondo problem. (See also Refs. jBipSj .) 

C. Strong dot-dot coupling case 

Next we discuss electron transport when Vc/A > 1. The solid line in Fig. [|(a) represents V g dependence of G; the 
solid and dotted lines in Fig. |](b) represent V g dependence of K and A, respectively (Vc/A = 10. 0, J/A = 1.0). For 
this case, the gate voltage dependence of G is fairly robust against J except that the region of G = for low gate 
voltages. When G — 0, k — J/2 and A = (spin-singlet state in the dots) as in the case of Vc/A < 1. Both of k 
and A change gradually against the gate voltage (there is a simple relation of 2k/ J = 1/2 — A/A, see Appendix C). 
These results contrast to the ones for the weak dot-dot coupling case. 

The robustness of G against J is explained as follows. At the top of the broad peak of G, n ~ 1 as seen in the 
previous section. On the right-hand side of the peak, n < 1 and spin-spin coupling J is not relevant. On the left-hand 
side of the peak, n > 1 and the antiferromagnetic interaction becomes more effective with decreasing V g . However, the 
transition from the Kondo states to the spin-singlet state takes place at quite a low gate voltage because the Kondo 
temperature 7k [Eq. (|l6|)1 is much larger than 7k, i for strong dot-dot coupling. In addition, t/A is always larger 
than unity when Vc/A > 1. Thus the spin coupling J does not introduce qualitatively new features in expression 

©■ 

D. Intermediate region of Vc/A 

Finally, we show V g dependence of G by changing the magnitude of Vc/A in Fig. |[ As Vc increases, (i) the height 
of the plateau increases, (ii) the sharp peak of G broadens, and (iii) the position of the sharp peak of G shifts to larger 
V g . When Vc/A > 1, the peak position is fixed at V g ~ 0. At the peak of G, n = 2 when Vc/A <§; 1 while n = 1 
when Vc/A 1; n decreases gradually from 2 to 1 as Vc increases. In conclusion, electron transport and electronic 
states in the dots change continuously from the weak to strong dot-dot coupling regions. 

V. ASYMMETRIC DOUBLE QUANTUM DOTS 

In this section, we discuss the effects of the energy difference AE = El — Er between the dots. First we describe 
the case of Vc/A < 1. In Fig. M, the gate voltage dependence of G is shown for AE = (thin dotted line) and 
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AE (thick dotted line) in the absence of J. These plots show that AE has no effect on the height of G when 
the gate voltage is low enough and both of the dots accommodate one electron. This is because the dot-lead coupling 
makes the Kondo resonant state at the Fermi level in the leads, irrespective of the position of the dot levels. The 
single peak of the transmission probability T{lo) is located at the Fermi level, as discussed in Fig. ||(a). As a result, 
G is insensitive to AE. In the presence of J, the gate voltage dependence of G is shown for AE = (thin solid line) 
and AE ^ (thick solid line). The sharp peak of G appears, but the peak height is less than 2e 2 /h. 

Next, we describe the case of Vc / A > 1. In Fig. [?], the gate voltage dependence of G is shown for AE = (thin 
dotted line) and AE ^ (thick dotted line) in the absence of J. For low gate voltages, til = ur = 1 and G is 
independent of AE as in the case of Vc/ 'A < 1. In the valence fluctuating regime, where n decreases as V g increases, 
ni, 7^ tir, and hence, electron transport depends on AE prominently. The broad peak of G splits into two peaks. 
They are less than 2e 2 /h in height. At the left-hand peak, decreases from unity to zero with increasing V g while 
rifj ; ~ 1. At the right-hand peak, ur decreases from unity to zero while til — 0. These peaks are understood by 
the Kondo resonance between a molecular orbital in the dimer and leads. One molecular orbital has the amplitude 
mainly in the left dot, and the other in the right dot. The peak heights of G are determined by the amplitude of 
the orbitals in the other dot. Note that the gate voltage dependence of G is fairly robust against J as in the case of 
AE = (compare the thick dotted and thick solid lines in Fig. ^). 

Finally, we show the intermediate region of Vc/A in Fig. [|. When Vc / A <C 1, G has a sharp peak. As Vc/A 
increases, the sharp peak broadens and a weak peak develops at V g ~ AE/2. The width of the weak peak is ~ A. 
The line shapes of these two peaks become similar to each other. Electron transport changes continuously from the 
weak to strong dot-dot coupling regions, as discussed in Fig. H in the case of AE = 0. 



VI. CONCLUSIONS AND DISCUSSION 



We have investigated the conductance through coupled quantum dots in series, as a function of the gate voltage, 
below the Kondo temperature. Electron transport is characterized by the ratio of Vc/A. When Vc/A < 1, a plateau 
of G appears because each dot forms the Kondo resonant state with conduction electrons in a lead and G is determined 
by the hopping between the Kondo states. For low gate voltages, they are replaced by the spin-singlet state and G = 0. 
A sharp peak of G appears with the height of 2e 2 /h, which is due to a coherent electron transport brought by the 
cooperation between the Kondo and spin-singlet states. When the energy levels are different in the dots, the sharp 
peak survives, but the peak height is less than 2e 2 /h. When Vc/A > 1, the molecular levels between the Kondo 
resonant states appear. The gate voltage dependence of G has a broad peak. The antiferromagnetic spin coupling is 
ineffective, compared with the case of Vc/A < 1, because (i) the Kondo temperature of a coupled dot is large and (ii) 
the transmission probability always has two peaks. When the energy levels are different in the dots, the conductance 
G has two peaks, reflecting asymmetric molecular levels. 

We should mention the validity of the mean-field approximation. First, the mean-field theory predicts the first-order 
transition between the Kondo and spin-singlet states when Vc/A < l/7rp3 it is an actifact of the mean field theory, 
because the NRG results show that when Vc ^ 0, the transition is always smooth.cJ : cJ For the present problem of 
coupled dot systems, this artifact appears on the left-hand side of the sharp peak of G; a discontinuous jump of G is 
seen for Vc/A = 0.2 in Fig. Second, in the valence fluctuating regime, the fluctuations around the mean field of 
the slave boson are not negligible. Hence, our calculations are not enough for quantitative discussion. 

In real coupled quantum dots, the intersite Coulomb energy, Ui n tULnR, is present. Then the ground-state energy 
E g . int is given by 

E sM t = E ss + C/ int (l - 6|)(1 - b 2 R ), (20) 

in the mean-field approximation. Equation ( po| ) indicates that the intersite Coulomb energy raises the dot levels El 
and En effectively. As a result, the valence fluctuating regime is extended and the Kondo regime shifts to lower V g . 
Note that we have assumed that the on-site Coulomb interaction U is infinite, and thus, {7; n t <C U . 

When [Tint is large enough, .the coupled dots cannot accommodate more than one electron. This situation has been 
investigated by Pohjola et al!Ei. In this case, the coupled dots is regarded as a single quantum dot with two levels 
when the dot-dot tunneling is large. Pohjola et al. have found a splitting of the Kondo resonance due to the level 
difference AE between the two dots. The one peak is on the Fermi level Ep and the other is located at Ep + AE. 
This is in contrast to the present case in which the split peaks are located below and above Ep. 

Finally, we comment on an asymmetric dot-lead coupling case, ^ A^. This asymmetry introduces the following 
two effects. First, the maximum of the conductance is less than 2e 2 /h by a geometrical asymmetry 4:AlAr/(Al + 
A/j) 2 .0 Second, the difference of dot-lead coupling <5A = (A^ — Ar)/2 introduces a difference in the number of electrons 
between the dots, resulting in El ^ Er. Electronic states are, therefore, asymmetric; the situation is similar to the 
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case of El ^ En discussed in Sec. |v|. [Equation (AS) in Appendix A is modified by replacing A — > A + b\ SA and 
Ai — * Ai + 6§ <5A. The latter modification means that ul ^ ur when SA ^ 0.] 
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APPENDIX A: DERIVATION OF THE GROUND STATE ENERGY 

BUS 



First, we review a single impurity model with the impurity level E 
Hamiltonian 7i 



H="£E(k)cl m c km +J2Eflf n 



Vb 



We start from the following mean-held 

E (4m/™ + f^km) + A(6 2 - 1) (AI) 



with E = E + A and b is the mean field value of the slave boson. 
The free energy F of the system is given by 



F = -| J duj ln(l + e- 0u )p(u) + \{b 2 - 1), 



(A2) 



where p(uj) is the density of states, which is given by the retarded Green function G r (to) of electrons in the dot 

1 



/?H = — G r (u>) = — 3 — . 

7T 7T (jj - E - E(w) 



(A3) 



with the self-energy T,(u) 



h 2 V 2 



1 b 2 V 2 

— = —in n(uj) = -iA$>(w). 

N *-r< w - E(k) +ir] N v ' 



(A4) 



Here p(u>) is the density of states of conduction electrons in a lead, p = p(uj = 0), $(w) = p(u>)/p, and A = %pb V /N 



To derive Eq. (A4), we have as sum ed th at p (uj) varies smo othly when compared with A 



By the substitution of Eqs. (A3) and (A4) into Eq. (A2) and integration by parts, the free energy (A2) becomes 



N 

F = — Im 

7T 



duf{w)$(u}) In (w - f) + A(& 2 - I), 



(A5) 



with £ = E + iA. 

If the leads has a wide flat band of width 2D, $(cj) = 1(— D < w < D), the free energy is 



N f 

F = — Im / dw/(w) In (w - + X(b 2 - I). 

7T J-D 



At zero temperature, this yields 



F = ^ImULU)-l)\ + M!P-\l 



An extension to the coupled dot system is straightforward. The free energy is 



F=-JmJ2 / dw/(a;)ln(w-ep)+ V K{b 2 a ~ 1 



JV 



(A6) 



(A7) 



(A8) 



a— L,R 
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with 



£ P = E + iA + PJ (^i + jA i) 2 + f2 



(A9) 



Here B ,i = (Pl ± E R )/2, A ,i = &o.iA witri & o,i = ( & l ± & !?)/ 2 > and t = v cVK ~ b f + K - £p are tne eigenvalues of 
the matrix given by 



V c + K E R + iA R 



Then the ground-state energy E gs = jVe gs [= F(T = 0)] is written as 

e gs = 2(P - E*)(b 2 - 1/2) + 2(Ei - AE) b\ + -k 2 + - V Im 



P=± 



(A10) 



(All) 



where E* = V g + A/7rln {Dn / A) with the bandwidth D in the leadsB In experiments, E* can be controlled by the 
gate voltage V g (we redefine E* as V g ) 



Minimizing the ground-state energy (All), we determine five parameters, b a , A Q (a = L, R), and k, self-consistently. 
The self-consistent equations are 




iA + 



PlltVc 



In 



7r£i 



= 0, 



2 Eh. 



AE\ 1 



bt-b\^(E 1 +iA l Y+t^ 
iE 1 A-bl(A* + V* + KVcl<fbt-b\) ni ^ 



Pin- 



(P x + iAi) 2 + f 2 



0. 



2 (& 2 -l/2) + i^Im 



In 



0. 



26? + -^Im 



(Pi+iAi) p]n *& 
(Ei+iA^+t* A 



2k 1 



(Pr+iA^ + t 2 



=Pln 



7rff 



= 0, 



0. 



(A12) 

(A13) 
(A14) 
(A15) 

(A16) 



APPENDIX B: EXPRESSION OF CURRENT 

In this appendix, we derive the expression of the current using the real-time functional integral method. The 
alternative derivation is given in Ref. ^5|, which is based on the equation of motions for the Green functions. We take 
h = 1 in this appendix. 

All of the physical quantities can be obtained in principle by calculating the time evolution of the density matrix 
Pit) 

p(t) = U\t)p init U{t), (Bl) 
where pi n a. is the initial density matrix and U (t) is the time evolution kernel of the Hamiltonian H 

U(t)=Texp(-iH(t)), (B2) 

using the time ordering p rod uct T . 



Corresponding to Eq. (Bl), we introduce the generating functional W[j 1 , j2],3EI 

exp (W [j'x, j 2 ]) = Tr Ht)p init U j *(t), (B3) 
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where U^{t) is the time evolution kernel of the Hamiltonian H' J = H — j-f — f-j with j = (jLm, jam), f = (fLm, fn.m) 
and their conjugates j and /. Equation (B3) is-,e*a,luated using the real-time functional integral method, which is 
the functional integral on the Keldysh contour;c3~Ell the suffix "1" and "2" in j correspond to the upper and lower 
branches of the Keldysh contour, respectively. 

For the Hamiltomaii (|l0|), Eq. (B3) is given by the following functional integral on the degrees of freedom of 
electrons in the dotsH 



exp(W[J]) = / VfVf exp ( ^ / duj F G" 1 F - J ■ F - F ■ J 

,m=±l/2 ' 



(B4) 



F = (/ + ;Lm(w), f-;Lm(u), f+;Rm(u), f--R m (w)) and J = (j—,Lm{w), j+;LmM, j_;R m (cj), j+;R m (u)) with /+ = 

(fi + f- = fi — /21 j+ = (ji + j2)/2, j- = ji — 32- The matrix G" 1 is the inverse of the Green functions: 



G- 1 H = 



( G-t(co) t \ 

G-{(u;) A L f L {u) t 

t G£r(w) 

V t G-^uj) AnfR^J 



(B5) 



with the Fermi function f a (u>) of the electrodes with the chemical potential [i a and the retarded and advanced Green 
function G~\ a (a — L, R) of electrons in the dots: 



G-i(w) = uj + E a + iA ai 
G-}Ju) = oj + E a -iA a . 



(B6) 
(B7) 



After integrating over the degrees of freedom of electron in the dots, the generating functional W is finally given by 
G(oj), 



W = - Yj / diJ 3 G H J 

m=±l/2 

The current 7 ([ll]) is rewritten as 

I = -ieV G ^ / du ((/ +; Lm/+;flm(w)) - (f+;Rmf+;Lm{u))) 



This is represented using the generating functional W 

5 5 



I = —ieVc / duj 



^J--L m 5J--R m 6J- ; R m 8J-.Lt 

= -2ieV c J duo [Gi,3(w) - G 8 ,i(w)] , 
where G n>m {ui) is the (n,m) component of G. 7 is finally given by 

7 = 2e / dw T(w) - 



exp(iy[J+, J_ 



(B8) 

(B9) 

(BIO) 
(BH) 

(B12) 



with T(u>) in Eq. 



To derive Eq. (B13), we have assumed A Q , E a , and k are determined by the self-consistent equations for the 
thermal equilibrium case \il — (J-r- Thus, the present calculation is valid only in the limit of ul — * [J-r- Extensions 
to the case of finite V s d is possible within the framework of the present method. See Refs. |3(i| , [43| . 
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APPENDIX C: KONDO TEMPERATURE WHEN Vb/A » 1 



by 



In this appendix, we derive the gate voltage dependence of A when J = 0. The self-consistent equations are given 



f^-i/2) + ^£ Im 



In 



TT^p 

A 



= 0, 



(PU c +*A)ln( ^ 



(CI) 
(C2) 



where £ P = E + PV C + iA with Vc = b\V c and A = b\A. 



When V g < —Vc/2, the Kondo state is formed in each dot and lead and thus E = 0. Then Eq. (C2) is reduced to 

V c 



A 

V g arctan(Vc/A) H In 



0. 



(C3) 



resulting in the first line of Eq. (0). 

When — Vc/2 < V g < Vc/2, the broad peak of G appears. For this case, the bonding level is located near the 
Fermi level in the leads and the antibonding level is far away from the Fermi level. Then we can use the following 



two conditions: (i) E = Vc ( Re£_ = 0) and (ii) Yip I m 
substituting Eq. del into Eq. (Jc|) yields 



In 



A 



Pin 



A 



= Im 



In 



A 



2V C 



2tt 



J] Re 



In 



17tA 



In 



7T 2Vb + iA) 



Then 



(C4) 



The logarithmic terms are negligible because n = 1 so that A ~ A/2 and Vc ~ Vc/2. Thus the gate voltage 
dependence of A reduces to the second line in Eq. (p"7|). 

When Vc/2 < V g where n < 1, the coupled dot system is essentially the same as the noninteracting system; n is 
determined by the position of the energy level: 



n= 1 



arctan (V g / A) 



(C5) 



This result yields the third line of Eq. (|17|), because A = A(l - n/2)/2. 

Fin ally, we mention the relation between n and A when J / (Sec. M). The self-consistent equation for n, Eq. 



2k 1 v^ t 
— + -> Im 



Pin 



A 



= 0, 



is simplified for —Vc/2 <V g < Vc/2. Using the conditions (i) and (ii), we obtain 

2k A 



(C6) 



(C7) 
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FIG. 1. A quantum dot dimer connected in series. 



FIG. 2. The transmission probability T(u>) as a function of (u) — E)/A with the energy of incident electrons ui, the level 
position E and width A of the Kondo resonance. El = Er and J = 0. (a) Vc/A = 0.3, (b) 1.0, and (c) 10.0. Arrows indicate 
the positions of the Fermi level in the leads for several values of the gate voltage V g /A. 

FIG. 3. (a) The conductance G versus the gate voltage when Vc/A = 0.3. El = Er. The dotted and solid lines represent 
the cases of J = and J/ A = 9.0 x 10~ 4 (U/A = 4 x 10 2 ), respectively, (b) The Kondo resonance width A (dotted line) and 
the spin-singlet order parameter k (solid line) versus the gate voltage for J/ A = 9.0 x 10 -4 . The inset shows G (broken line), 
A (dotted line), and k (solid line) as functions of the gate voltage, near the sharp peak of G. 

FIG. 4. (a) The conductance G versus the gate voltage when Vc/A = 10.0. El = Er. The dotted and solid lines represent 
the cases of J = and J/A = 1.0 (U/A = 4 x 10 2 ), respectively, (b) The Kondo resonance width A (dotted line) and the 
spin-singlet order parameter k (solid line) versus the gate voltage for J/A = 1.0. 

FIG. 5. The conductance G versus the gate voltage for several values of Vc/A from 0.2 to 3.0. El = Er. All the curves are 
vertically offset for clarity. 



FIG. 6. The conductance G versus the gate voltage when Vc/A = 0.3. For J = 0, AE/A = 0.0 (thin dotted line) and 2.0 
(thick dotted line). For J/A = 9.0 x 10~ 4 (U/A = 4 x 10 2 ), AE/A = 0.0 (thin solid line) and 2.0 (thick solid line). 



FIG. 7. The conductance G versus the gate voltage when V c /A = 10.0. For J = 0, AE/A = 0.0 (thin dotted line) and 8.0 
(thick dotted line). For J/A = 1.0 (U/A = 4 x 10 2 ), AE/A = 0.0 (thin solid line) and 8.0(thick solid line). 

FIG. 8. The conductance G versus the gate voltage for several values of Vc/A from 0.2 to 3.0. AE/A — 2.0. All the curves 
are vertically offset for clarity. 
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